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Abstract 

The azimuthal anisotropy of high-p^ particle production in non-central heavy ion colli- 
sions is among the most promising observables of partonic energy loss in an azimuthally non- 
symmetric volume of quark-gluon plasma. We discuss the implications of nuclear geometry 
for the models of partonic energy loss in the context of recent RHIC data and consequences 
for observation of jet quenching at the LHC. 



^Contribution to tiie CERN Yellow Report on Hard Probes in Heavy Ion Collisions at the LHC. 



1 Introduction 



In order to interpret data on nuclear collisions from current experiments at the Relativistic Heavy 

Ion Collider (RHIC) and future experiments at the Large Hadron Collider (LHC), it is necessary 
to have knowledge of the initial conditions. There are large uncertainties in the estimates of the 
initial produced gluon density, Pg{To) ~ 15 — 50/fm^ in central Au + Au at y/s = 130, 200 AGeV 
and Pg{TQ) ~ 100 — 400/fm^ in central Pb + Pb reactions at ^/s = 5500 AGeV, since widely 
different models (e.g. see [1, 2]) seem to be roughly consistent with data [3]. It is, therefore, 
essential to check the energy dependence of the density of the produced quark-gluon plasma 
(QGP) with observables complementary to the particle multiplicity dN'^^/dy and transverse 
energy dEx/dy per unit rapidity. High-^^T observables are ideally suited for this task because 
they provide an estimate [4] of the energy loss, Ai?, of fast partons, resulting from medium 
induced non-Abclian radiation along their path, as first discussed in [5, 6] in the context of 
relativistic heavy ion reactions. The approximate linear dependence of Ai? on pg is the key 
that enables high-p^ observables to convey information about the initial conditions. However, 
/S.E also depends non-linearly on the size, L, of the medium [7, 8] and therefore differential 
observables which have well controlled geometric dependences arc also highly desirable. 

A new way to probe A.E in variable geometries was recently proposed in Refs. [9, 10]. 
The idea is to exploit the spatial azimuthal asymmetry of non-central nuclear collisions. The 
dependence of /S.E on the path length L{(f)) naturally results in a pattern of azimuthal asymmetry 
of high-pr hadrons which can be measured via the differential elliptic flow parameter (second 
Fourier coefficient), V2{pt) [H, 12]. Before we show the sensitivity of the high-py V2ipT > 2 GeV) 
to different initial conditions we briefly discuss the various model calculations for the "elliptic 
flow" coefficient V2- 

1. The elliptic flow parameter V2 was first introduced in the context of relativistic hydrody- 
namics [12] and reflects the fact that due to the macroscopic sizes of large nuclei many 
aspects oi A-\- A collisions are driven by nuclear geometry. In non-central collisions the 
interaction region has a characteristic "almond-shaped" form as shown in Fig. 1. 



Hydro calculations convert the ellipticity of the reaction volume into momentum space 
azimuthal asymmetry 

, _ - jy') ^ ipl) - (pD ^ = /o^^^^°^^^^¥iSw^ ... 

through the higher pressure gradient along the small axis. The elliptic flow is thus perfectly 
correlated to the reaction plane and can be used for its determination [13]. Hydrodynamic 
simulations [14, 15] typically describe well data from relativistic nucleus- nucleus collisions 
at v^AfAT = 200 GeV up to ~ 1.5 — 2 GeV and it is not unlikely that at LHC energies of 
V^NN = 5.5 TeV the region of validity of those calculations may extend to — 5 GeV. 

2. Initial conditions can also be mapped onto flnal state observable distributions by solving 
covariant Boltzmann transport equations as in cascade models (partonic, hadronic, and 
multi-phase). Elliptic flow in this approach is generated via multiple elastic scatterings. 
Calculations are sensitive to the choice of initial conditions [16] and are currently limited by 
statistics to ~ 6 GeV. It is interesting to note that they can match the high-pT behaviour 
of the V2 but require extremely large initial gluon rapidity densities dN^ /dy ~ 16000 [16] 
and/or string melting [17]. 

3. Memory of the initial parton density, reaction geometry, and the consequent dynamical 
evolution is also retained by large transverse momentum partons (and fragmented hadrons) 
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through their jet quenching pattern [9, 10]. While this approach is discussed in more detail 
below, it is important to emphasize here that at the single inclusive jet (or hadron) level 
the resulting high-pr azimuthal asymmetry is also perfectly coupled to the reaction plane. 
It has been suggested that in the limit of very large energy loss the momentum asymmetry 
is driven by jet production from the boundary of the interaction volume [18]. 

4. Recently, a classical computation of the elliptic flow at transverse momenta > in 
the framework of gluon saturation models has been performed [19]. It was found that 
the azimuthal asymmetry is generated already at proper time r = 0, i.e. it is built in 
the coherent initial conditions. The resulting elliptic flow coefficient was found to vanish 
quickly V2ikT) oc k^'^{R-'^ - R'^) above (~ 1 GeV for RHIC and ~ 1.4 - 2 for LHC 
energies) which is not supported by the current data. 

5. An approach that does not associate azimuthal asymmetry with the reaction plane has also 
been presented [20]. Both high-pT and low-pT V2 emerge as a back-to-back jet correlation 
bias (with arbitrary direction relative to the reaction geometry for every pT bin) . For large 
transverse momenta V2 oc lnp7^//i suggest an easily detectable factor of 3 increase in going 
from pt = 5 GeV to pT = 100 GeV at LHC. The pT-integrated V2 oc 1/Qs at LHC exhibits 
~ 50% reduction relative to RHIC. (It can also be deduced that V2 is larger at the SPS in 
comparison to RHIC in this model.) 

The methods for V2 analysis can be broadly divided in two categories: two-particle meth- 
ods discussed, e.g., in [21] and multi-particle methods [22, 23]. In two-particle methods the 
error on the determined V2 from non-flow (non-geometric) correlations is 0{l/{v2M)), where 
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M is the measured multiphcity. With multi-particle methods this error goes down typically to 
0(1/(^2-/^^)), i.e., smaller by a factor of order M. Although it is not possible to completely 
eliminate the non-flow components to V2, experimental techniques based on higher oder cumu- 
lant analysis [22, 23] will be able in many cases to clearly distinguish between between reaction 
geometry generated azimuthal asymmetry and back-to-back jet bias. 



2 Energy loss in a longitudinally expanding plasma 

We first generalize the finite energy gluon bremsstrahlung theory [8, 24, 25, 26] to take into 
account the expansion (neglected in [9]) of the produced gluon-dominated plasma while retain- 
ing kinematic constraints important for intermediate jet energies. The GLV reaction operator 
formalism [26] expands the radiative energy loss formally in powers of the mean number, 
of interactions that the jet-|-gluon system suffer along their path of propagation through dense 
matter. For a jet produced at point xq, at time tq, in an expanding and possibly azimuthally 
asymmetric gluon plasma of density p{x,t), the opacity in direction v{(f)) is 

POD 

X{<t>)= dT a{T)p{xo + v{(fy{T -to),t) . (2) 

J TO 

Note that the gluon-gluon elastic cross section, (7(r) = Q'K0i^g/2n'^^f{T), and the density may 
vary along the path. The explicit closed form expression for the n**^ order opacity expansion 
of the gluon radiation double differential distribution for a static medium is given in [26]. The 
generalization to the case of arbitrary medium similar to Eq. (2) is trivial and amounts to 
replacing the static weight with the nested position integrals in the case of a dynamically evolving 
quark-gluon plasma: 
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Fortunately, the opacity expansion converges very rapidly, and the first order term was found to 
give the dominant contribution. Higher order corrections decrease rapidly with the jet energy 
E. All numerical results discussed here include 2'^'^ and 3'''* order correlations between scattering 
centers that can alter the leading order dependence outlined analytically below. 

In the case of 1-l-lD Bjorken longitudinal expansion with initial plasma density po = p('R)) 
and formation time tq, i.e. 

/^W = Po(^)", (4) 

it is possible to obtain a closed form analytic formula [10] (under the strong no kinematic 
hounds assumption) for the energy loss due to the dominant first order term [8]. For a hard jet 
penetrating the quark-gluon plasma dAE^^^/Edx oc f {Z {x , t)) dr / \{t) , where x ~ lj/E is 
the momentum fraction of the radiated gluon and the formation physics function f[Z[x,T)) is 
defined in [10] to be 

/(^, r) = r [1 - COS ( uZ{x, r))] . (5) 

^0 u{l+u) 

With Z{x,t) = (r — to)ii'^{t) /2xE being the local formation physics parameter, two simple 

2a. 

analytic limits of Eq. (5) can be obtained. For x ^ Xc = p{to)'^Tq^' L^~~ /2E = L//^(L)/2£J, 
in which case the formation length is large compared to the size of the medium, the small 
Z{x,t) limit applies leading to f{Z) ~ ttZ/2. The interference pattern along the gluon path 
becomes important and accounts for the the non-trivial dependence of the energy loss on L. 
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When X <^ Xc, i.e. the formation length is small compared to the plasma thickness, one gets 
f{Z) InZ. In the x ':§> Xc limit [10] the radiative spectrum reads: 

dx ~ 2(2 - a) A(ro) x ' ^ ' 

The mean energy loss (to first order in x) is given by 

The logarithmic enhancement with energy comes from the Xc < x < 1 region [26]. In the case 
of sufficiently large jet energies {E oo) this term dominates. For parton energies < 20 GeV, 
however, corrections to this leading Inl/xc expression that can be exactly evaluated numerically 
from the GLV expression and are found to be comparable in size. The effective AE/E in this 
energy range was found to have a plateau [27]. We note that for Bjorken expansion with a = 1, 
the asymptotic energy loss can be expressed in terms of the initial gluon rapidity density as 

It is evident from Eq. (8) that the appropriate variable to drive non-Abelian jet energy loss 
calculations with is the effective initial gluon rapidity density dN^ /dy. 

3 Parton energy loss and nuclear geometry 

For nucleus-nucleus collisions the co-moving plasma produced in an A + B reaction at im- 
pact parameter b at formation time tq has a transverse coordinate distribution at mid-rapidity 
Pg{Y-,z = 0,ro). In studying jet production and propagation in nuclear environment it is not 
always technically possible to perform the Monte-Carlo averaging over the jet production points 
coincidentally with the simulation of parton fragmentation. It is therefore useful to separate the 
medium dependence of the mean jet energy loss as a function of the extent of the nuclear matter 
traversed and the azimuthal angle (f) relative to the reaction plane. In the Bjorken, Eq. (4), 
and linear f{Z) « 7rZ/2, Eq. (5), approximations the total energy loss is proportional to a line 
integral along the jet trajectory r(T, (p) = r + v{(j)){T — tq), averaged over the distribution of the 
jet production points 



T/i(r) = / p/i (r, z) and TyiB (6) = / d^rr^(r)rB(r — b) depend on the geometry. In particular, 
for a sharp uniform cylinder of radius J?efF one readily gets r^(r) = (A/7ri?gg)^(i?eff — |r|) and 
Tab{0) = A'^/'^RcS- therefore define the effective radius of the sharp cylinder equivalent 

to a diffuse Wood-Saxon geometry via F{0, (/))wood-Saxon = -^(0, ^^')sharp cylinder- For Au + Au 
collisions and a = 1 the above constraint gives R^s ~ 6 fm. 

For a non- vanishing impact parameter h and jet direction v{(t>), we calculate the energy loss 

as 

^E{hA) Fib,4>) AEjO) _ AEjO) 

= m^) = ^^^'"^^ ' ^^^^ 

where the modulation function R{b, (/)) captures in the linearized approximation the b and (f) 
dependence of the jet energy loss and also provides a rough estimate of the maximum ellipticity 



4 



DC 



1.4 
1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



A+A (A~200) 
Pb+Pb, Au+Au 



Wood-Saxon geom. 
Sharp cylinder geom. 



b = 2fm 




0.0 



2.0 4.0 

(|) [rad] 



6.0 



Figure 2: The modulation function R{b, (j)) is plotted versus </> for several impact parameters b = 
2, 6, 10 fm from Ref. [10]. Diffuse Wood-Saxon versus uniform sharp cylinder geometries are compared. 
The most drastic difference between these geometries occurs at high impact parameters. 



generated via correlations to the reaction plane. Fig. 2 shows the R{b,4>) modulation factor 

plotted against the azimuthal angle (j) for impact parameters & = 2, 6, 10 fm. Note that R{h, (p) 
reflects not only the dimensions of the characteristic "almond-shaped" cross section of the inter- 
action volume but also the rapidly decreasing initial plasma density as a function of the impact 
parameter. 

In order to compare to data at pT < 2 GeV at RHIC and pT < 5 GeV at LHC, one must 
also take into account the soft non-perturbative component that cannot be computed with the 
eikonal jet quenching formalism. The hydrodynamic elliptic flow [12] was found in [14] to have 
the monotonically growing form V2s{pt) ~ tanh(pT/(10 ± 2 GeV)) at ^/s = 200 AGeV and 
to be less sensitive to the initial conditions than the high-p^ jet quenching studied here. The 
interpolation between the low-pT relativistic hydrodynamics region and the high.-pT pQCD- 
computable region can be evaluated as in [10]. 



Fig. 3 shows the predicted pattern of higli-p^ anisotropy. Note the difference between sharp 
cylinder and diffuse Wood-Saxon geometries at 6 = 7 fm approximating roughly 20-30% central 
events. While the central (6 = 0) inclusive quenching is insensitive to the density profile, 
non-central events clearly exhibit large sensitivity to the actual distribution. We conclude that 
V2{pt > 2 GeV, b) provides essential complementary information about the geometry and impact 
parameter dependence of the initial conditions in A + A. In particular, the rate at which the 
V2 coefficient decreases at high pT is an indicator of the diffuseness of that geometry. Minimum 
bias STAR data at RHIC [28, 29] for > 6 GeV now seem to support the predicted [9, 10] 
slow decrease of V2 at large transverse momenta. Recently in [30] hadron suppression in Au + 
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Figure 3: The interpolation of V2{pt) between the soft hydro dynamic [14] and hard pQCD regimes is 
shown for 6 = 7 fm adapted from Ref. [10]. SoUd (dashed) curves correspond to sharp cylindrical (diffuse 
Wood-Saxon) geometries presented in Fig. 2. 



Au {Pb + Ph) relative to the binary scaled p + p result at — 5 GeV for RHIC conditions 
iV^NN^dN^ /dy) was found to be approximately equal to the quenching factor at LHC at a 
much larger transverse momentum scale pT — 50 GeV. One may thus anticipate proportionally 
large (~ 10 — 15%) azimuthal asymmetry for high pT at the LHC. 

4 Energy loss of jets in transversely expanding medium 

Transverse expansion can only be very crudely modelled in Eq. (4) by taking a > 1. Such 
approach neglects the time needed for the rarefaction wave to propagate from the surface to the 
center of the interaction region. To derive an improved analytic expression taking transverse 
flow into account, we consider next an asymmetric expanding sharp elliptic density profile the 
surface of which is defined by 

+7.^^ = 1. (11) 



(i?^ + ■u^r)2 {By + VyTf 

The area of this elliptic transverse profile increases with time, r, as A±^{t) = Tr{Rx + VxT){Ry + 
VyT). The plasma density seen by a jet in direction (pQ starting from xq = (xo,yo) inside the 
ellipse with a specified initial condition tqPo = l/{TrRxRy) dN^/dy is 

1 dN^ fl 

p{r,(f)o;xo,yo) ■■ 



IT dy \T J \Rx + VxT J \Ry + VyT 
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\ {R. + V,Tr ^ {Ry + VyTf ) ■ ^ ' 

We approximate the assumed 4>q independent screening /x(r) k, gT{T) = 2(/9(r)/2)^/^ since 5 c± 2 
and p = (16C(3)/7r2)r3 ~ 2T^ for gluon plas ma and define r(^o) as the escape time to reach 
the expanding elhptic surface from an initial point xq in the azimuthal direction (pQ. We take 
a;(^o) = 2t(0o) (p(''"(?!'o))/2)^^^ to estimate an upper bound on the logarithmic enhancement 
factor. A short calculation leads to: 

4 RxRy dy ttx - ay cl)(0o) 

where ax = Vx/Rx,ay = Vy/Ry. This expression is a central result for transversely expanding 
media and provides a simple analytic generalization that interpolates between pure Bjorken 
1+lD expansion for small ax^yT, and 3+lD expansion at large Oj; „r. 

In the special case of pure Bjorken (longitudinal) expansion with Vx = Vy = 

In this case, the energy loss depends linearly on r(^) and we recover the result of Eq. (8). 
Another special case is azimuthally isotropic expansion with Rx = Ry = R and Vx = Vy = vt- 
Taking also the longitudinal Bjorken expansion into account leads in this case to 



We note that for a jet originating near the center of the medium and fully penetrating the plasma 
the enhanced escape time due to expansion t = R/ {1—vt) compensates for the 1/(1+vtt(^o)/-R) 
dilution factor. Therefore, in this isotropic case, the extra dilution due to transverse expansion 
has in fact no effect of the total energy loss: 

AE[^^{b = fm) ^ AE^^^{b = hn) , (16) 

modulo logarithmic factors which become sizable only for large vt- An important consequence 
is that the inclusive azimuthally averaged jet quenching pattern in central collisions is approxi- 
mately independent of transverse expansion. We have checked numerically that Eq. (16) holds 
for realistic transverse density profiles [14]. 

In non-central collisions, the azimuthal asymmetry of the mean energy loss can be expanded 
in a Fourier series and characterized as 

A£:J^(0) = AE{l + 2S2{E)cos2(P + ---) . (17) 

It is correlated to the final measured elliptic "flow" of jets and hadrons and has been evaluated 
by using a full hydrodynamic calculation from Ref. [14]. In this case we use the parameterization 
eBC of [14] to initialize the system and treat gluon number as conserved current to calculate the 
density evolution needed in the line integral Eq. (9), where it replaces the naive Bjorken (ro/r)" 
expansion. We average over the jet formation points the density of which is given by the number 
of binary collisions per unit area as in the Woods-Saxon geometry used in Ref. [10]. We find that 
the azimuthal asymmetry of the energy loss is strongly reduced for realistic hydrodynamic flow 
velocities. This implies a much smaller V2 at high pT than obtained in Ref. [10] where transverse 
expansion was not considered and poses questions about the observability of the effect at LHC. 
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5 LHC-specific remarks 

There are several important aspects in which LHC and RHIC wih differ significantly. We briefiy 
discuss the implications of those differences for high-pj^ V2 measurements: 

1. Currently at RHIC at ^/s = 200 AGeV the > 2 — 3 GeV regime is perturbatively 
computable [30] (modulo uncertainties in the baryonic sector [31]). At LHC the px region 
which is not accessible through the pQCD approach may extend to transverse momenta as 
high as 5-10 GeV. This would imply the validity of the relativistic hydrodynamics in this 
domain, the extent of which can be tested by looking for marked deviations in the growth 
of V2{pt), saturation, and turnover. 

2. Estimates of the initial gluon rapidity density at LHC vary from dN^ /dy = 2500 to 
dN^/dy = 8000. This would imply very large parton energy loss, at least in some regions 
of phase space. In this case jet production for moderate transverse momenta may be 
limited to a small shell on the surface of the interaction region, leading to a constant 
V2{pt) purely determined by geometry [18]. 

3. Since mean transverse expansion velocities at RHIC have been estimated to be on the 
order of vt — 0.5 through relativistic hydrodynamics fits, it is natural to expect even 
larger values at LHC. This may lead to a significant reduction of the observed azimuthal 
asymmetry as discussed above. An important prediction of the approach put forth in [10] 
is that V2{pt) exhibits a slow decrease with increasing transverse momentum. This can be 
used to distinguish azimuthal anisotropy generated through energy loss from alternative 
mechanisms. 



5.1 Jet impact parameter dependence at the LHC 



In light of the discussion in Sec. 5 it is important to asses the feasibility of azimuthal asymmetry 
measurements for largc-E'T jets via detailed simulations. The impact parameter dependence 
of jet rates in Pb + Pb collisions at the LHC was analyzed in [32]. The initial jet spectra at 
■\/s = 5.5 TeV were generated with PYTHIA5.7 [33]. The initial distribution of jet pairs over 
impact parameter 6 of ^ + ^ collisions (without collective nuclear effects) was obtained by 
multiplying the corresponding nucleon-nucleon jet cross section, crj^^, by the number of binary 
nucleon-nucleon sub-collisions [34]: 



d^ 



^ib,v~s) = TAAiby^Vv~s) 



(18) 



with the total inelastic non-diffractive nucleon-nucleon cross section cr™^ ~ 60 mb. 

The rescattering and energy loss of jets in a gluon-dominated plasma, created initially in 
the nuclear overlap zone in Pb + Pb collisions at different impact parameters, were simulated. 
For details of this model one can refer to [32, 35]. To be specific, we treated the medium as a 
boost-invariant longitudinally expanding fluid according to Bjorken's solution [36] and used the 
initial conditions expected for central Pb + Pb collisions at LHC [37, 38, 39]: formation time 
To 0.1 fm/c, initial temperature Tq ~ 1 GeV, gluon plasma density pg « 1.95r^. For our 
calculations we have used the collisional part of the energy loss and the differential scattering 
cross section from [32]; the energy spectrum of coherent medium- induced gluon radiation was 
estimated using the BDMS formalism [7]. 

The impact parameter dependences of the initial energy density eo a-i^d the averaged over 
(fi jet escape time {tl) from the dense zone are shown in Fig. 4 [32]. (tl) goes down almost 
linearly with increasing impact parameter b. On the other hand, eo is very weakly dependent 
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of b {Ssq/eq < 10%) up to b on the order of nucleus radius Ra, and decreases rapidly only 
at 6 > Ra- This suggests that for impact parameters b < Ra, where ~ 60% of jet pairs are 
produced, the difference in rescattering intensity and energy loss is determined mainly by the 
different path lengths rather than the initial energy density. 

Figure 5 shows dijet rates in different impact parameter bins for E';^* > 100 GeV and the 
pseudorapidity acceptance of central part of the CMS calorimeters, jr/J'^*! < 2.5, for three cases: 
(i) without energy loss, (ii) with collisional loss only, (iii) with collisional and radiative loss. The 
total impact parameter integrated rates are normalized to the expected number of Pb+Pb events 
during a two week LHC run, i? = 1.2 x 10^ s, assuming luminosity L = 5 x 10^^ cm~^s~^. The 
maximum and mean values of dN'^^^'^^ / db distribution get shifted towards the larger 6, because 
jet quenching is much stronger in central collisions than in peripheral ones. Since the coherent 
Landau-Pomeranchuk-Migdal radiation induces a strong dependence of the radiative energy loss 
of a jet on the angular cone size [40, 41], the corresponding result for jets with non-zero cone size 
^0 is expected to be somewhere between [Hi) {9q 0) and {ii) cases. Thus the observation of a 
dramatic change in the 6-dependence of dijet rates in heavy ion collisions as compared to what 
is expected from independent nucleon-nucleon reaction pattern, would indicate the existence of 
medium-induced parton rescattering. 

Of course, such kind of measurements require the adequate determination of impact param- 
eter in nuclear collision with high enough accuracy. It has been shown in [42] that for the CMS 
experiment the very forward pseudorapidity region 3 < [ryl < 5 can provide a measurement of 
impact parameter via the energy flow in the very forward (HF) CMS calorimeters with resolu- 
tion ai) ~ 0.5 fm for central and semi-central Pb + Pb collisions (see details in the section on jet 
detection at CMS). 

5.2 Jet azimuthal anisotropy at the LHC 

While at RHIC the hadron azimuthal asymmetry at high-pT is being analyzed, at LHC energies 
one can hope to observe similar effects for the hadronic jet itself [35] due to the large inclusive 
cross section for hard jet production on a scale Et 100 GeV. 

The anisotropy of the energy loss (AE) goes up with increasing b, because the azimuthal 
asymmetry of the interaction volume gets stronger. However, the absolute value of the energy 
loss goes down with increasing b due to the reduced path length L (and Eq at b > Ra, see 
Fig. 4). The non-uniform dependence of AE on the azimuthal angle ip is then mapped onto 
the jet spectra in semi-central collisions. Figure 6 from [35] shows the distribution of jets over 
ip for the cases with collisional and radiative loss (a) and collisional loss only (b) for 6 = 0, 6 
and 10 fm. The same conditions and kinematical acceptance as in Fig. 5 were fulfilled. The 
distributions are normalized by the distributions of jets as a function of ip in Pb + Pb collisions 
without energy loss. The azimuthal anisotropy becomes stronger in going from central to semi- 
central reactions, but the absolute suppression factor is reduced with increasing b. For jets with 
finite cone size one can expect the intermediate result between cases (a) and (b), because, as we 
have mentioned before, radiative loss dominates at relatively small angular sizes of the jet cone 
^o(— > 0), while the relative contribution of collisional loss grows with increasing 6*0. 

In non-central collisions the jet distribution over (f is approximated well by the form A{1 + 
Scos2(^) , where A = 0.5(iVmax + A^min) and B = (iVmax - A''mm)/(A''max + ^^min) = 2 (cos2(^). 
In the model [35] the coefficient of jet azimuthal anisotropy, = (cos 299^°*) event , increases 
almost linearly with the impact parameter b and becomes maximum at 6 ~ 1.2Ra- After that 
v^2^ drops rapidly with increasing b: this is the domain of impact parameter values where the 
effect of decreasing energy loss due to the reduction of the effective transverse size of the dense 
zone and the initial energy density of the medium is crucial and cannot be compensated by 
the stronger volume ellipticity. Anther important feature is that the jet azimuthal anisotropy 
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Figure 4: The impact parameter dependence of the initial energy density en{b)/eo{b = 0) in nuclear 
overlap zone (solid curve), and the average proper jet escape time (tl) /Ra of from the dense matter 
(dashed curve) [32]. 
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Figure 5: The jct+jct rates for Ej^* > 100 GoV and jryj"*! < 2.5 in different impact parameter bins: 
without energy loss (solid curve), with collisional loss (dashed curve), with collisional and radiative loss 
(dotted curve) [32]. 
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Figure 6: The jet distribution over azimuthal angle for the cases with collisional and radiative loss (a) 
and collisional loss only (b), > 100 GeV and jryj"*! < 2.5 [35]. The histograms (from bottom to top) 
correspond to the impact parameter values 6 = 0, 6 and 10 fm. 



decreases with increasing jet energy, because the energy dependence of medium-induced loss is 
rather weak (absent in the BDMS formalism and ~ In in the GLV formalism for the radiative 
part at high Et)- 

The advantage of azimuthal jet observables is that one needs to reconstruct only the direction 
of the jet, not its total energy. It can be done with high accuracy, while reconstruction of the jet 
energy is more ambiguous. However, analysis of jet production as a function of the azimuthal 
angle requires event-by-event measurement of the angular orientation of the reaction plane. The 
methods summarized in Ref. [11, 43, 44] present ways for reaction plane determination. They 
are applicable for studying anisotropic particle flow in current heavy ion dedicated experiments 
at the SPS and RHIC, and may be also used at the LHC [35]. Recently a method for measuring 
jet azimuthal anisotropy coefficients without cvent-by-event reconstruction of the reaction plane 
was proposed [45]. This technique is based on the correlations between the azimuthal position of 
jet axis and the angles of hadrons not incorporated in the jet. The method has been generalized 
by taking as weights the particle momenta or the energy deposition in the calorimetric sectors. 
It was shown that the accuracy of the method improves with increasing multiplicity and particle 
(energy) flow azimuthal anisotropy, and is practically independent of the absolute values of 
azimuthal anisotropy of the jet itself. 
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6 Conclusions 



The azimuthal anisotropy of high-pT hadron production in non-central heavy ion colhsions is 

shown to provide a vahiable experimental tool for studying both gluon bremsstrahlung in non- 
Abelian media and the properties of the reaction volume such as its size, shape, initial parton 
(number and energy) rapidity densities, and their subsequent dynamical evolution. The satura- 
tion and the gradual decrease at large transverse momentum of the reaction geometry generated 
V2, predicted as a signature complementary to jet quenching of strong radiative energy loss in a 
dense QCD plasma [10], seem now supported by preliminary data extending up to pT — 10 GeV 
at RHIC. 

The initial gluon densities in Pb + Pb reactions at ^/sj\!pf = 5.5 TeV at the Large Hadron 
Collider are expected to be significantly higher than at RHIC, implying even stronger partonic 
energy loss. This may result in interesting novel features of jet quenching, such as modification 
of the jet distribution over impact parameter [32] in addition to the azimuthal anisotropy of the 
jet spectrum. The predicted large cross section for hard jet production on a scale of Et ~ 100 
GeV will allow for a systematic study of the differential nuclear geometry related aspects of jet 
physics at the LHC. 
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